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I. INTRODUCTION 


A. BACKGROUND AND RELATED WORK 


The study of electromagnetic energy transmission is but 
one important area in microwave and millimeter-wave engi- 
neering, where the electromagnetic waves are travelling 
through some transmission medium, which provides’ the link 
between the transmitting and receiving part of a transmis- 
Slon system. 

In recent years, fin-line has gained in importance as a 
transmission medium in millimeter wave circuit constructions 
[Rem 1 - 4]. Fin-line has been found superior to micro- 
strip at millimeter wavelengths as the former provides eased 
production tolerances, low dispersion, broad single mode 
bandwidth, moderate attenuation, better compatibility with 
hybrid devices, greater freedom from radiation and higher 
mode propagation, comeinedewrtem the ability to construct 
Simple transitions to conventional rectangular waveguide. 

Figure 1.1 shows a 3-dimensional view of fin-line. The 
Sstucture may be viewed as a slotline with a shield, a ridged 
waveguide with dielectric or a slab loaded waveguide with 
Finse 

The fin-line structure was first proposed for millimeter 
wave integrated circuits in 1974 by Meier [Ref. 1]. An 
early paper by Meier described the propagation mode as a 
variation of dominant mode in ridged waveguide. His proce- 
dure requires a test meaSurement to determine the equivalent 
Cioleetrie constant of the fin-line structure. This is both 
expensive and time consuming. Knorr and Kuchler [Ref. 5] 
presented a frequency dependent hybrid-mode analysis of slot 


line with open boundary using the spectral domain technique 





Figure 1.1 3-Dimensional View of Fin-line Structure. 
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which was suggested by Itoh and Mittra [Ref. 6], in 1975. 
Subsequently, the paper of prime importance in the estab- 
lishment of spectral domain technique for analyzing the 
fin-line structure, so-called Shielded slotline, was 
pesented by Knorr and Shayda in 1980 [Ref. 3]. One of the 
advantages of the spectral domain approach is that it is 
numerically more efficient than the conventional methods 
that work directly in the space domain. This is due prima- 
rily to the fact that the process of Fourier transformation 
of the coupled integral equations in the space domain yields 
a pair of algebraic equations in the spectral domain that 
are relatively easier to handle. Another important advantage 
is that the Green's function takes a much simpler’ form in 
the transform domain, as compared to the Space domain where 
no convenient form of the Green's function is known to 


exist. 


Boe PURPOSE 


In solving the open boundary slotline problem, only 
exponential functions arise and there are no numerical prob- 
lems during computation. For a closed boundary structure, 
however, hyperbolic functions are required. Knorr and Shayda 
[Ref. 3] found that overflow and underflow problems resulted 
during numerical computations. To eliminate these problems, 
extensive algebraic manipulation of the spectral domain 
equations is required. 

This thesis presents the method to calculate the charac- 
teristic impedance of fin-line without overflow and under- 
flow in equations. Therefore, this thesis is a direct 
SxEentromuon Knorr s work. This work is accomplished by lots 
of algebraic manipulations. In particular, since hyperbolic 
Sine and cosine functions cause overflow and underflow 


errors, the equations developed by Knorr need to be put in a 


ih 


form where only the hyperbolic tangent function appears. The 
characteristic impedance 1S computed after the spectral 
domain technique to find the dispersion characteristic. 
Following the theoretical analysis, an explanation of 
the computer program used in determination of characteristic 
impedance is presented. Numerical results are then compared 
with known data for ridged waveguide, slab loaded waveguide, 
slot line and fin-line [Ref. 3]. This comparison estab- 


lishes the accuracy of the numerical results. 


eZ 


IIT. THEORETICAL ANALYSIS OF FIN LINE 


A. FIELD AND BOUNDARY CONDITIONS 


The fin-line supports a hybrid 


nants of TM and TE modes are then 


Beka pO) a 


Hz 


kz gh (x, 4) ef 


where the scalar potential 


Helmholz equation, and we 
that [7 = +)j8. 
Further 
a A ae 
Ke, = Ki — 6 
zs 2 ‘ 
with ky = WMAER , We eee gS 


finline as stated in Figure 1.1. 


Through Maxwell's curl equations 


field. 


functions 


assume lossless 


The axial compo- 


(eqn 2.1) 
(eqn 2.2) 


the 


ZS dn Seine tans 


propagation so 


(eqn 2.3) 


defining spatial region of 


the transverse field 


components are then determined by these axial components and 


can be given as 


e h 
= (72 - jou) el 
ey (PF + yond) el 
He= (P +jweES a 
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(eqn 2.4) 


Ceci) 


(eqn 2. oc) 


z Pa = eqn : 
Hy= CP a - we 3) e™ iieniiiin 


where propagation in the z direction is assumed. We will 
also assume here that €,= &% = €& and €)= €€r. 

At y = h, + D 
Applying boundary conditions at the walls in region 1, 


tangential field components must be zero. 


Ex (YX, h+D,Z) =O (eqn 2.8) 
Fx\ CF ht+D,z) =O (eqn 2.9) 
At y = D 
At the interface between region 1 and region 2, tangential 


field components must be continuous. 


Ez (%,D,2) = Ezz GSE, a) (eqn 22, 10) 


Ex (X, 0,2) = Exz CX, 0, 2) (eqn 2.11) 


Also the electric fields at y = D will exist only in the 
Sot 


0 Ixjz 2 
Ezy X,D,z) = 
rz Ww. 
C2) € IMh< a (eqn 2.12) 
° xj 2 
Ex, (4, 0,2) = WW) 
exnel™ Ae ooa (eqn 2.13) 
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Similarly, tangential magnetic fields must be discontinuous 


by corresponding surface current densities. 


% 


oe xyz 
Hay (10,2) Hea (K,002) = } me 


O Ixic & 
(eqn 2.14) 
; rz 
Me™ |x z¥ 
Hx (A, 0,2) ~ Hx2.(K,0,2)=] 9 a 
Oo Xl < = 
(eqn 2.15) 


AT y = 0 
The tangential field components at the interface between 


region 2 and 3 must also be continuous. 


E22 (40,2) = E23 (K%,°,2) (eqn 2.16) 
Ex2 (K,0,%) = Ex3 (4, 0, 2) (ccnmeaion) 
H22 (K, 0, 2) = H23 (%, 0) 2) (eqn 2.18) 
Hxo (X,0,2) = x3 (K, O%) (eqn 2.19) 
mee y =e tlgo : 
Once again at the shield wall in region 3, the tangential 


field components must be zero. 


F23CX,-h2,2) =O (eqn 2.20) 


Ex3 CX hat ee) =e (eqn 2.21) 


IES 


eae / 2 
The final boundary conditions occur at x = b/2 where the 


tangential components must be Zero in all regions. 


Ee: (tbb,4,4) =0 (ean 2.22) 


Exar (+/2, 4,Z) =O (eqn 2.23) 


B. SPECTRAL DOMAIN APPROACH TO DISPERSION CHARACTERISTC 


The scalar potential functions can be transtormed ante 


the spectral domain via Fourier transform defined as 


{blah = baaw= {i bauer” dy, 


(eqn 2.24) 


The scalar potential functions satisfy Helmholz's equations 


in the three spatial regions, thus 
2 ae 
Ve oi ale Kei p; = 0 (eqn 2.25) 


where denotes the two dimensional Laplacian operator in 
the transverse (x,y) direction. The Helmholtz equation 2.25 


are eranstormed irmtoe 


3; be¥) = (ola — Kes) $: Gm, 4) 


34 (eame2. 26) 
ole 2 a: we 2 
where yz = lm - Key = olm + © - K . Above equation has solu- 
tions after applying boundary conditions at y = D+ hy, and 


yo = = Tig 


B, dm, 4) a AS en) stmh y\ COT (eqn 2.27) 
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$5 Gn, 4) = B° Ga) stmh Ya + C&En)cash Ys} (eqn 2.28) 
5 (on , ¥) = D& Gm) stmh y; Chat 4 ) (eqn 2.29) 


BM Gm) 4) = Ah dn cosh ¥; (D+h-Y) (eqn 2.30) 


$2 (in, 4) = B" Gn) sinh A + ch@n)cosh ¥Y (eqn 2.31) 
Bz mn, 4) = Dh (lm cosh Y3Chat J ) , (eqn 2.32) 

ee 
ame ob" evem (eqn 2.33) 


Xn = 


enn ph Bad 


(eqn 2.34) 
The following observation about the solutions in region 
2 must be emphasized. Any wave on this inhomogeneous wave- 
purges Sermeture will partly travel through air and partly 
through the dielectric slab. It is important to observe at 
this point that Yj; may be less than zero in any of the 
three regions of the structure under certain conditions. 
When m= 0 and ky approaches kg (where kg is the wave 
number for free space), 6 is less than kx; and so tae O. 
Under this condition the hyperbolic functions in all three 
regions are replaced by trigonometric functions. If Ko < 6 < 
kz, then y,; and Ya Premeteater thanezero and Y; < 0 for some 
values of and the trigonometric functions replace _ the 
hyperbolic functions in the spatial region 2 only. This 
Suggests that the nature of the field is dependent upon the 
values of the transform variable, we EQ the, conditions 


2 


when Y: ew Y; replaces y; such that (yan = -)y,. The 


- 


wv; 


eight unknown coefficients A through D are not indepen- 


dent, but can be related to each other through the conti- 


nuity conditions of the field components at the interfaces 


between the three spatial regions. 
If we denote the Fourier transforms of x- and z-directed 


current density and electric field component by 


Ex Gin) = FL exw} (eqn 2.35) 
Ean =F {ea} (eqn 2.36) 
Jy Ee) aaa {5(00| (eqn 2.37) 


jaca) =F {20} ; (eqn 2.38) 


The resulting set of linear equations may be written in 


matrix formas follows: 


AS Oo 
Be O 
(Me) 4) 2s 
ch 
p’ Ee (econ 
Ac ° 
Be 0 
[Melee | oy 
A 4 
p . 
(eqn 2.40) 
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The matrices [ME] and [MJ) differ in only the last two rows. 


Each is a square 8 x 8 matrix. Using equations (2.39) and 
(2.40), we may write 
@) © 
Oo oO 
_| = ° 
(MsJ[Me] | - | = | - 
ww ~ 
Ex Sx 
~~ fw 
Ez Jo/. 
(eqn 2.41) 
From equation (2.41), using the four elements in the lower 
right hand corner of the matrix MJ ME , we obtain 


GE m,P) Fe n,e)] { Ex@m)]  ( TeG) 


Gs, dm, 6) Ga a,B) E2 Gn) Jz Ga) - (eqn 2.42) 


wolere= Chiiee2 xX 2 matrix [c} SQaealismeune FOUrler Eranstorms 
of the components of the dyadic Green's function for this 
SVnUucruLe. 

A solution to equation (2.42) is obtained using the 
Method of Moments [Ref. 7]. For this problem, we have 


chosen to approximate the field between the fins as shown in 


Bteure 2.1; 


w 
Xl o> 
Ex = 
O else where Ceah 2s) 
C2 (K) -oOo. (eqn 2.44) 


eo 


ed 


-w/2 +w/2 


Figure 2.1 Assumed Electric Field Component in Slot 
in x-direction versus x for Fin-line. 
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The dispersion problem is now reduced to the form 


EF Gab) [ExAe)|” = 0 


(eqn 2.45) 
where 
AL ae Joln XK simCnW/2) 
E dn) = ‘, el) ET dK = AW 


Pemunertcdlm@isearch for the value of which satisfies equa- 
tion (2.45) yields the propagation constant for the dominant 
fin-line mode. 


From the equation (2.39) 


my A© + MBE + mi3sc® =O 


(eqn 2.46) 
MaAS + Ma2 Bo + mo3c& + 
nas AD + Mt BN + Man Cl = O (eqn 2.47) 
ons ne = pres (ean 243) 
mn hee wines AM = DE (eqn 2.49) 
m33C& + m34 DS = O Canes) 


maz Ce + mada Do + m4 B+ mag D" = 0 (Cones 
psn ch + ms De = 0 (ep 2. 5) 
Mer BE + me.De + M67 Ch 4 sD = 0 (eqn 2.53) 


Haommeciations (2.50) and (2.52) 


Ze 


3 4- 


he - TE (eqn 2.55) 


Using equations (2.48) and (2.49) 


“= (aE) 


(eqn 2.56) 


Wee \ aw aon eB 
A = (Ont ) OE - ae : (eqn, 2.78) 


Substituting in equations (2251!) and@@7e5> eve scbpEeam 











Go h 
os ene + ae ae ¢ (eqn 2 sen 
ph = (MEE? a yc + (= ae. 
(M34. (eqn 2.599 
Substituting equations (2.48), (2.49), (2.52), a (2.53 \eee 
equations (2.46) and (2.47) we obtain 
ea Q\2 Mi es 
(DI ie AW ior mE) Bre « 








Ma} 
-Cor) (aa mf) DE 


(eqn 2.60) 
c= (oi) ne SE + Gor) ) (eae 
ae) B Dé + (Er) (Ge) Ee (ect) ca) 
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where 


[D| = 4 Q22 - Audi2 


= j Cduda2 - dardiz) 20) som hd) cosh( 20) - 


We define 





Ou = [ mis + yz ( ea 2 )| 


(Mb2. M34 


=~ dy cosh (%,P) 


pie ) — “er 


= (rma ( 
Gia ( - M62 M58 M™ 62. 


= ) di2 (LP) simh (Y2D) 


Q2\ = (ea lf 


62 M34. 


M44 M33 M43 
Mob ag M34 | MAL ) 


= da cosh ( YP) 


dar = [ may + tag ( Mas Me) 


£ Mss 
M3 Moy AMET 
Riis (M62. MCs M62 ) 
= j daz CaP) smhOhD) 


Alsou we detine 


Ze 


(eqn 


(eqn 


(eqn 


(eqn 


LO.) 


pai) 


2.64) 


Daisy), 


; (WED) OBOY Ckeaoy 
dy = - (Ke2d) [\ BA (WED) OLD) Cke3P)* 


(7,0) tanh CLD) | 
C¥3D) tanh [C¥30)( 449) 


7 Gnd) (6D) 
din = (Kad) [oe ee a4 = 


(eqn 2.66) 


(eqn 2.67) 


a (WED) (3D) ¢ Kad) 
day = — Gn0)(80) C (WE2D) (LD) ——_ 
CYP) 4amh (2D) 
(Yb) tamh (0s) (p)] —_— (eqn. 2 69) 


2 (Ker D130) tom COST) 
doa = (upoy[ C+ SAC Od donb OLD) 


(Amd) CPPS” 
tr Soe C= ] ; 


(WO MAD) (WEP) C¥2DYz (eqn 2.69) 


Nomalized constants are defined as follows 


(Kei) = (Keo = 20Y 01 -CYKF] CAY 
(Kady = (ax [Er — (WNP) COA 
WMD = 240 T* (P/n) 


WED = WE3D = Zz CP/n) 
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Weap= Er/bo CD/A) 
Bp = 20 OA)Q/X) 


M 20 2 oh even 
~ ce (2) bh odd 


UY = CVF = God + (AF CCAS 1} CAS 


Op = adi + OUP COWS ErJ CA). 


C. CHARACTERISTIC IMPEDANCE 


The definition of the characteristic impedance for an 
ideal TEM transmission line is uniquely given by static 
quantities. Since the fin-line supports a hybrid mode, no 
unique definition of the characteristic impedance can be 
found. 


A useful definition, however, is 


ue 
2Fav4 (eqn 2.70) 





Zo = 


where V, is the slot voltage defined as 


hy 
ae \ oy AdAx =|. (eqn 2.71) 


22) 


e, (x) is arbitrarily selected as 1/W so that W*ey(x) = 1. 


Pavg 1s given by 
Pavs ae + Re { EXH* - dz da 


= 4 Re \\ (Exthy® = Ey Ht) andy 


Geqn 2.72) 


Pasval's theorem is applied to eq(2.88) to obtain 


pe 


Oth 
Pavg = REG Zaye Jy, CEx Hen YS Gn, 4D 


This expression must be evaluated in each of the three 
regions of the fin-line shown in Figure 1.1 Therefore the 


power flow may be expressed by 
Pav = 35 = CPitPa + Psd (eqn 2.73) 


Since Pav4 can be determined after finding the value of 
KIN Equations including only hyperbolic tangent functions 
and the slot voltage Vg have been developed the character- 
istic impedance. The lengthy algebraic manipulations are 


shown in appendix B. 
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III. COMPUTER PROGRAMMING 


A. NUMERICAL ANALYSIS 


The 
upon the 


computation of characteristic 


solution to the dispersion 
for the fin-line under consideration. 
wave propagation constant 
known before any other investigations can be 
only in this case are the 
transform domain known. 

ength ratio A/a is already developed by Prof. 
of the 


for the numerical evaluation of the time average 


next task is the preparation 


For ease in numerical calculations and for 


purposes all geometric parameters are normalized 


impedance is 


In other words, 


or wavelength ratio X/a 


Knorr. 


based 


chracteristic problem 


the 


must be 


Started since 
scalar potential functions in the 


The computer program for the wavel- 


The 


appropriate equations 


power flow. 
programming 


as follows; 


h;y/D ; fin location relative to the positive "y" side 
wall normalized with respect to D 

h2/D ; fin location relative to the negative "y" side 
wall normalized with respect to D 

b/D ; waveguide height normalized with respect to D. 


ie 


infinite integration in 


is observed that for the power 
equation (2.88) 
infinite Series due to the 
variable km. 
coefficients of theSe Series 
with respect to the variable 
be saved by Summing over a half interval only. 
The time average power flow equations are 
six different cases as follows; 
case l ; (CY, Din < 0 in region l 


case 2 ; (Y, DY > Q in region l 


flow computations 


Preliminary numerical investigations 


the 


is replaced by an 


discrete nature of the transform 


of the 


indicated an even distribution 


so that computation time can 


prepared for 


case 3 ; Gyo < 0 in region 2 and ( Y3D - < Q 


Zy 


in region 3 


case 4 ; opis < 0 “an regionyZseand ( yD y > 50 
in region 3 

case 5; (y, DY~ > 0 in region 2 and ( Ys D y < QO 
in region 3 

case 6 ; (¥2D)~ > O in region 2 and (Y; D ig > O 


in region 3 

After this preparation of the equations for the computer 
programming, the only remaining task. is to investigate the 
numerical . integration with regard to its limits and the 
procedure to be used. One of the limits of the summation 
discussed previously extends to infinity so it iS necessary 
to determine an appropriate at which to truncate the 
computations. 

From the representative examples shown in Figure 3.1, it 
is seen that the coefficients in this infinite series are 


found to decay rapidly so that a finite approximation yields 


good results. Figure 3.1 shows the characteristic impedance 
as a function of n, the number of terms in the truncated 
series, for various slot widths. The value n= 50 is 


sufficiently large to be considered infinite for all prac- 
tical purposes. It is noted that when w/b = 1, Ex(dm) = 0 
for n > QO. Therefore, whenever w/b = 1 only the n= O term is 


computed in impedance calculations. 


B. COMPUTER PROGRAMMING AND LIMITATION 


The first part of the FINIMP computer program finds the 
f 
A/A which makes equation (2.45) equal to zero. The assump- 
tion is that A/A is between 0.1 and 3 in the giga hertz 


range and continuously decreases as_ the frequecy is 
increased. Since this program is very sensitive there is a 
possibility that the wrong value is searched. In that case, 


we can easily notice because there is abruptly jumped value. 
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Zo 


At that time, if the assumed range of X/n is resticted to a 
smaller interval than before the correct value can be 
obtained. 

The rest of the FINIMP computer program is the computa- 
tion of characteristic impedance. Since all equations are 
composed of hyperbolic tangent functions, there iS no over- 
flow in the computer program. But if (Y4D)~ < QO tangent 
functions replace the hyperbolic tangent functions. At that 
time, there is a possibity of overflow because tan nf/2 is 
infinite fomwne= ele wee . In that case the tangent 
function should be set to some value which is the maximum 


value of computer ability WG (1-16) > | taniGe ile 
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IV. NUMERICAL RESULTS AND COMPARISONS 


To check the accuracy of the numerical results gererated 


by the computer program, comparisons are made with data 
available in the literature for ridged waveguide, slab 
loaded waveguide, slotline, empty waveguide and fin-line as 


outlined in [Ref. 3]. 


A. RIDGED WAVEGUIDE 


When = 1 and w/D < 1 or when the dielectric substrate 
thickness Dis reduced to zero, the fin-line structure 
becomes ridged waveguide with zero thickness ridges. 

The impedance of the ridged waveguide has the following 
relations [Ref. 8], 





Zz es Zo00 
o> 

(1 - a ybule ; (eqn!) 
When the guide width andthe slot width are equal, the 


ridged waveguide becomes an ordinary rectangular waveguide 


for which 


a Mike 2n 
om NN Ev& ~~ ob ia (eqn 4.2) 





When the guide width is increased at fixed slot width, the 
characteristic impedance at infinite frequency and the free 


Space cutoff wavelength go asymptotically to 
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Characteristic Impedance Z vs. 


Frequency for a Ridged Waveguide. 


Figure 4.2 


S3 


Ze 2b 
oe EyEo A Zama (eqn 4.3) 


The cutoff wavelength A, was determined from [Ref. 11] for 
several values of W/b. 

Equation (4.1) was then used to calculate the character- 
istic impedances and these values were compared with the 
results using spectral domain method of the computer program 
discussed in chapter], . Two results are very close 
together. Figure 4.1 and 4.2 shows the computer results. 
These figures agree with the previous Knorr and Shayda's 


CPesuisese 


B. DIELECTRIC SLAB LOADED WAVEGUIDE 


A variation of the fin-line structure of Figure 1 where 


w/b = land €y > 1 results ina dielectrie slao loaded 
rectangular waveguide. The slab loaded waveguide has also 
been studied by Vartaniam et  aleiRet wo They consider a 


guide with the slab centered in the waveguide and they 
present an analytical expression for the voltage impedance 
Zpy at the center of the slab. The impedance computed in 
this analysis was specified at the edge of the dielectric 
slab. Therefore a relationship between the two impedances 
must be defined before any comparison can be made. The impe- 
dance at the edge of the slab can be easily defined by the 


expression 
ae 


a Fx d~ +) 
or Pv F center 
Xx (eqn 4.4) 
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SLAB LOADED WAVEGUIDE 
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where Ege is the field at the edge of the slab and Bee 
is the field at the center of slab, thus 


a= Zoey cos (2) 


Cece) 


where the various quantities are defined in [Ref. 9] as 


Se C/2 (eqn 4.6) 





(Ly = Ey (Ay- a) (eqn 4.7) 


a (20)°(F aes = Can earme)) 


and c = Dis the dielectric slab thickness. 

The characteristic impedance of a slab loaded guide is 
computed uSing the spectral domain method. The results were 
compared with equation (4.5) results and [Ref. 3]. Good 
agreements are obtained. Figure 4.3 shows the’ slab loaded 


wave guide impedance 


Gre SLOTLINE 


If w/D < 2 and Ey is sufficiently high for the fin-line 
Structure of Figure 1.1, the field is tightly bound to the 
Slot. For this condition the presence of the shield will 
have little effect if the walls are sufficiently far removed 
from the slot. In this case the fin-line will behave like a 


Slotline. This behavior is illustrated in Figure 4.4 where 


a, 


the characteristic impedance of a fin-line with w/D=1, &, 
= 20 have been plotted. Also these results are in good 


agreements with [Ref. 3]. 


D. SEVERAL FIN LINE IMPEDANCE CURVES 


Fin-lines are generally enclosed with a shield that is 
compatible with the dimensions of the standard rectangular 
waveguides for the millimeter wavebands. Above 22GHZ, all 
these guides have aspect ratios b/a = 0.5. Further, the 
fins are most often centered in the guide and printed using 
D = 0.005 inch substrate with €y= 2.2. 

Figures 4.5 - 4.12 show the impedance computed for fin- 
line with WR(28), WR(19) and WR(12) rectangular waveguide 
Shields. These figures may be compared with data for the 
Same structures as presented in [Ref. 3]. Such a comparison 
Shows excellent agreement for all but the small values of 
W/b. For small values of W/b (W/b = 0.1, W/b = 0.05, W/b = 
0.02), the results of the FINIMP program show lower impe- 
dances than the results presented in [Ref. 3]. As the W/b 
is smaller, the difference is larger; for W/b = 0.1, 2 - 5 
ohm is lower, for W/b = 0.05, 9 - 15 ohm is lower and for 
W/b = 0702, 19 = 25 ohm iselowen: 
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V. CONCLUSIONS AND RECOMMENDATIONS 


A. CONCLUSIONS 


The spectral domain technique uSed in conjunction with 
Galerkin'’s method has been presented to calculate the char- 
acteristic impedance for the dominant mode of the fin-line. 
It has been shown that a matrix formulation of the problem 
permits the elements of the dyadic Green's function to be 
calculated. Solving the matrix equations'7 leads Eo 
containing only hyperbolic tangent function. These equa- 
tions circumvent the overflow and underflow problems which 
occurred using the previous formulation presented by Knorr 
and Shayda. 

Numerical results obtained using this method have been 
presented and compared to other existing data. Good agree- 
ment has been obtained in all cases thus’ establishing the 
accuracy and applicability of the method for the full range 
of structure parameters. 

There is a possibility that tangent functions may cause 
an overflow problem if (¥:D)~ is less than zero. In this 
case, the value of tangent function can be obtained within 
the capability of the available IBM 3033 used. 

In this thesis, paticular interest is devoted to the 
computation of fin-line impedance. Fin-line may exhibit the 
characteristics of ridged waveguide, dielectric slab loaded 
waveguide, slot lines, and conventional rectangular wave- 
guide. All of these structures are fin-line substructures. 
So, the computation of fin-line impedance permits all of 


these structures to be analyzed. 
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B. RECOMMENDATIONS 


Coupled fin-line will find future use in building direc- 
tional couplers and filters. For this purpose, the normal 
mode wavelengths and impedances are required. A special but 
important case is that of symmetrical lines. Coupled fin- 
lines for which the normal modes are odd and_ even. The 
program described here may be extended to cover this case by 
following the procedures outlined by Knorr and Kuchler for 
coupled slotlines [Ref. 5]. This should be accomplished to 


improve the utility of the program described in this thesis. 
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APPENDIX A 
SPECTRAL DOMAIN MATRICES 


The continuity conditions are transformed via equation 
(2.24) into the two dimesional Fourier domain. The solutions 
to the two Helmholtz equations piven by equations 
(meee eee) 06 6are§ «6Ssubseituted. Finally a matrix form of 


linear equations is derived as follows; 


Mi M2, O oO © O O AS oO 
May Maz Maz 0 mas MWLb Nth, O B oO 
Ceon Olas Misono Oo gO | lac O 
0 O m43 mM Oo Mag 0. AMBIT) DO] | O 
0 0 0 0 8 © sy Mss{| AR} |o 
O it o Mid © © Wey meg }}er}| |e 
mg OF OC oO tine on o Cc Ex 
Cc (OmmOuNTG) 0 anO LD Ez 
my Mo Ms © O O O © AS © 
M2) ‘tha Wa OO ws rhb m7 oO ge Oo 
O o 33 MR4 O ° O oO GC Oo 
Oo 0 mad 44 © Mus O MM4B11O |} | 0 
Pe et So) = 2 Se iiG 
one o m4 0 Oo my Mg 1] B? O 
oo 0 0 Ars A Mon o és i 
Ma, Mga. tgs 0 ngs my rig oO LD” Je 


The matrix elements of (Mg) and [(My} are normalized at 


this point with respect tcl. the dielectric substrate 
thickness. The normalized matrix elements are now presented 
in two forms. The first of the element equations is for 


De > O and the second is for (¥%D)> < 0. For the 
matrices [Me] and [My}, the elements my through mgg are the 


same. 


iol | CKerP) simh TO DCGY J 
| 
J Geof sim [eR’oy 69] 
ae — (Kad) simh (2D) 
i | - j(Kaoy sim 7/0) 
5 ~ Gedy cosh (LD) 
ad ~ (Kap} cos CYD) 
M4 = Mis = Me= Mi = Mg = O 
_ | (lm) (BP) Simh LODICMG)] 
J m0) (BO) sim Tev'odchypy J 
— ab) (BP) Simh Ch) 
~ J Gad) (6) sim CY") 
« ~— Gad) (BP) cosh CVD ) 
— @lnD) (¢P) cos (YD) 


| 


50 


J Comp) CYP) simh TONIC] 
— —3 (HMDIOY'D) simO N’p) ChYppJ 


f (wry YP) cosh C¥2P) 
Ma, = 
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(J (WME) CD) stm CBD) 
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_ - 7 (WEsD) C130) cosh COHPX HY) 


(W 
m (WE3 D)CYD) cos LOY DCD) 


mes = M™é6eb = O 

GlmD) 
nig = (OO 

Gn DCAD) 

- (lad) (6b) cosh COP] 
na | 

— Gnd) (pP) cos LOCH) 


me (dD) (BD) Simh CCN DIChYH)) 
J Gad) (BD) StmE CYDIChYHY] 
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m3, 2 en (72D) 
~}(Keaby sim CY¥/D) 
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34 
ii 
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APPENDIX B 
TIME AVERAGE POWER FLOW 


The following expression for the coefficients AS through 
DM are derived from equations (2.54) - (2.61). 
PemGy>)> < 0 
ot ot; 
> \ Capp sim LAV) 
qh = if (Ker DY DYEx — GnD) (80) D'Ez 
CWP) CKD CY’D) SIML YD hyp) 
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_ [Chet eo) 
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) det C.D) simh CD) CKeaD) cos (OP) (M46 


If (¥%,D) > 0 and (¥;D> 0 


60 


pf - [ WEsD)(aD Kea DY [4.08 - Dez 
GeDNOLD Teeb¥ femh LOCI detcosh(2D) | 


4. [Ga0(B0) [KeadY—(Kes DJ] Cage ee _ANDEx-dyo& 
CWE2D)OLD)CKe3 Dy ear a 


_ s [Gubd)@)[CKead)— ad) J) “de Ed DES 
a) (WUD)OSD) Kady  dekcoshQ2D) . 
+ C009 damh [ONC] If dud Ee “Ata 0 Es 
J (sp) (Kes eerie 


f= Kad) (dp 0 Ex —duD Es 
det cosh (Y2D) (Kea DY Smh (7730)C 4H) J 


he [ (Kea DY An 0 Ex — day D Ez) 
J L det QLD) smh QD) esp cosh [AX] 


From the above equations 
ct Cline < 0 


pe. _ Bs [1+ tow [DI H)] 


Dt (KeD)4 taney) Chip) ] 


Ahn KP Ex — WadXbDERT [ittan Lo’rxhoy) } 
Dt (WMO (DY Ker DF tan COOK YOY] 


Ae. (AY _ _ Ee {CK DY Ex — ad)(BD)Ez [i+ tan [alo] 


p* CWUDCN’D) (Kad * Aor? [0/06 J 
ary An : pe (Ay 
p* a 


[f(y oa 


ae.(aeyF ES Li- tanh TOD ChYp) I) 


Dt — (Key D)* tanh? (cx, orchypy] 


ah. (arit _ (Kadd& - aorGorEeT [i—tank coro) 
pt (WMD CY DY C10) Aaah LOND} 


pe. _ nuts DY Ex-GaP)(60) EaJLi- on LONGI] 
D+ JL CkeD+(W UP)CY,0) tank? Olin) Choy] 


acy. Ah uae AS. (AY 
prt pt 


Het OLD: < 0 


ec ae ey [ a ari (X09) 


ch {cr - [anes —dy Ea i “pte Ord). 2) 
det C’d) eS a. 
ce os (di2Ex—daoEs)(AnEx-du® passed itstant(%,b) 
J det caret 


Gee ch 7. ce. Gr 
p* 





ot 
en Gy AD 20 
= Ex _ jo 
= [- ey ; a *C1—dewnh (GD)] 
aac oer vex — \—stanh* CVD) 
e ak OD) i Saale) teal 3 


x wexK ~ z 
-~j(¢ LE, sete Ex -duy E: Osa ‘[- tank (yap) | ] 


CED) Aaah (BD) 





If (Y,D) < 0 and (¥3D)°< 0 


£ ot = (WE53D)( sD) (KeaDJ_ in ae 


Dt LW ED) ALO Kes 0F'C1;"0) tan(s0x'¥6)) 


(WED) OD) (Kes DY (¥3'0) ton (Cs DXHY] 


[ (dnd)(BD){CKe2D) ~ (Kes) 1] ; E ef) 
CWE2D)(¥2"P) (Kea 0) O* 


4 [Gadr(B0f Kead)= (esd9J i(2 | oF) 
(WED) (Y,/D)( Ka DY p* 


BN-(BAY _ Gav (BD)[K2D)—(Kead)] ‘{ EF ] 


e(é 
Dt (WMD) CD) (Ke3D)* pt 


_9 [ Glad (BD) (Ke2 BY — Kes D3 J , 
CWADICED) CKD 


64 


0) (Kead}tomCOs DICH) VEC ch 
C¥2’0) (Kes D~ II pt | 


. [ondesoTion Oso ehy) i CFT 


(1y’P)CKesd~ p* 


Bey i DOO) (KaD 
pt J LCWeD) (0) Keb ton (0D) CH] 


[ Gnd)(6D) (ker DJ —(Kes DY] cf. (eF 
(WUD)CND) (KexD Pa 


> [esd yO) (KeDY gions} {Ez ch : 


J we, 00h’) KesDYS LOD Kes0¥ JL j p+ 


~ [ GaD¥ (BDI [KsPI—CKes0JJ WP c& CH 
"J (WED) (Wd) Baur ll moe 





i [ GaP) (BD) chap) =a) 4 
JL WED OOK Kead* 


fas D) (Kea) tan(30) (6) ] ( ch. chyF 


CHD) CKD pt 


Be eer. (- (WES D\ 10)KeaDJ (<— ca 
(wes D)CL/D)(KesDY tan[077oX ra 
_ F. Gnd) GP [(KeaD) = Ga03 I oF) 
(WE2D)CVD)CKe3D)~ pt 


yk ch _ Be. 
p* D+ 


LOT CWE3D)(Y3"D) ken oa EP) 
pt esd) Kes0y-tan (0) 


— j FARGNCDY 919 ] (Fe 


(NED) (Y2'D)Ke3D~ 
Gre _ _ BF 
of Dt 


Ber os: eerie (ey 1 


(WMD) CY’D)( Kes D> 


66 


CY; 3 
+{ PgaTianll 06) ] s CF ) 


Cyc . pcr 


pt pF 








Br. (chy® ae GaDX ED KadY—(Kesd) 1 Tf c® chy 
pt (WMD) CY,’D) (KesD) j | 


3 3 
+; [9 een (3 YER ch. hyF ] 


cByech BAY 
-- 3 p+ 








Fee -(KeadDS ipieeettoin 6) c&. (LF 
_CKesb ton [(¥¥'0)%0)] { pt | 


pr. (oyF 
“= CADE i+ tout Cryo GE] 
(OF _ sas COED 1 eG 

D+ ea a som CC Y3' a [“ ) 0+ 


ry ph_ POY 
p* p* 


Lf (Y¥,D)- = 0 and seal iae = O 





_ (ne x 3 A 8 ; 
Be (Be _ pence 0)(¥30)(KesDY™ | E «ey 


p+ hid) (KesDytanh(CHD(MH) J L DA 


_9 [ (€30)C¥3D) CKea DY . 
WED)020) Kady tanh (Ospxhyp) J 


nD (G0) [Kadl- Kady Ji c& 
[ pa IL ce 


(dnd) Bd) [CKead)— Ces PV IP ch. (ch 
T CWED)CY0) (KesD It Was | 


Br. (phyr _ pabian lier (kesDY} \[ ee ] 
D+ CWUDICY/D Ks DY D* 


1 ; 2 GuD)(BD) L (Kea DI - (Kea) J . 
(Wak) (YD) (Ke3 D> 





ChockapT ath Lien) (re) 


+ [2 CeadF ooh LOCI) VP [ne cP) 
O10) (Ke3D)* 


ey (WEsD_)(¥sD)(Ke2 DY | 
pt (Wed) WD)( Kes DY tanh [(¥s0)(hp)I 


[ (dnd) BO) [(KeaDJ —CKes DJ] I c&. (<&F ; 
(WD) CY/D) Ke3D)* D* 


+ j [Wemaurckany fey (<9 (ch 
TC wep) 0) Kar Soca pt 


Ikea if ce. chy 
a) ea so : | 


Zs foe (Ckead~(KeaI} 7 
(WE2D)(1,’P) (KesD)* | 


(@ D)(Kaby tanh ((¥30%4)] [ay Cry ] 
C¥2’D) (Kea DY” 


69 


(ey ph Bt By 
pt p+ 


oars f (WEsD)(¥5D)(Ke2D ce. (ny 
a WE2 DI CAD) (KesD}* tom LODO VL j = 


_ -G.dx 60) (KesD) — (Kad). NE “ ( om 


CWEP)CYA'D) (KesD 


@yt.ch pe. 
po a 





be. (cer =o 
Dag J WE2D)( Y.D)(KesD tanh (C730) C*%)] 


“iF a0) BD) [(Kerd) (Kes DJ J (ec) 
J CWE2D) (YD) (keaD)* } * 


geyce pee 
pr pt 


BC _ pa 0) [CKead- (kes DI [et 
pt (Wud)CY'D)(KesD* ot 


—- GNsPEEO | (cet 


70 


(40)(KeaDY tanh [CYsDdC HI 17 (ceyF ch 
+ cee a ii ee ‘ 


(Bh). ceé 7 gh. (eye 
——_—— = pF 4 


car. PeaonponiTeant—Cesod (ct (ory 
pt = JL Cup ce’dX Kes jp4 


_ ~ FQD)Ke2d) tamh (040)(%) ] & , cr 
J Cyz0) (KesD¥ Dt 


pt pt 


zZ 


roy _ va | 1 ~ tanh’ C00) ('%)) [= 


D+ — LkspP JL tanh? [cy, ox) p+ 





ei [sey - sank cborton 


areal oes rset er 


71 


oop sp (OF 


p pH 


If (Y2D)* > O and (¥3D)~< 0 


Be (Be)™ _ (WEsDKYYD) (Kad) {jeer 
pt  LWEDICLD)(Ks0P tenn Lb pt J 


1 26 WEs0) (YSD) (Ker DY ] 
(WE2 PICLD)(Kead)” tam [(¥30)(4Y) 


GuD)(BO>[ (KeaDJ— Kes DF] [eh ch 
[ (WE2DICEDM Ke3Dy— D+ 


GnDD(BO)[Kead I — (Kes 03 J 77h (Af 
tL Cid) CEDED™ | : aT | 


BG _ dmD)( BD) El od 
, a (Wad) CYP) (Ke3P) p+ 


42 Gud) (BP) [Kady — e309) 7] 


CW AD)CI2D) (ke3D) 


ve epee [<r 


+ po eae De aot Repro \[ ch r | 


BEB _ = (esp vsyeKesdY 
p+ | Lr Caso hyp] 


— Kes }7 ce. OF 
mc 


- { (WesD) (VSD) (Ke2d)* (ey. ch 
TJ oh cel J D+ | 





_ + (lod) (80) [Kea - esd) I ce a 
YL WEDYWMDYCEDY Kad)t JL pt 


. CoP YON HDCP) 

til CWé2P)CYP)(Ke3DY> 
CVD) (Kea) kam [Ose ch. | 
CYP) Kez > Dt 


eer ph _ oe. UF 


p4+ Da 


Be YF Wes DY VHD) (KeadT- \< Cy 
pt (WEDICRD) (Kes sam(CysnHNJL D* | 


_. FdeDCB I-(KespI JP ch. ch 
i ae IE a ; 








B°.e [; (DEsDX VSD K Kad) [| 
D+ WEP) KKsPY tom (Xp UL D+ 


(KesD¥- (KeadT] 1 CeO* ch 
(WED) CYDXK3D J Dt 


(FY c&  pt.cyr 


p* D+ 
D+ 


(EY _ - CaP BP [Keds c&- (CF 
BY - it see er] 


74 


Ss ¥xD)C%) yee — 


(ey. me Da Bh (ec) 


pt p+ 


Bee y a = [eerGorl Great (KesPY] (Fe ch 
Coppy(py(kape JL spt 


- [RPT som (BOCA) ae oF 


CVD) (Kea DY 


p+  =——sOt 


“= 2 thon LOS hp CF 
. alk ert oe ee p+ 


=i) a | [i+ ta Lavor(tyoy) [EF 


pf. OF itdar LOD CHILE CY 
wat tom (CP) nee | lS 2) 


core pr _ of Or 
p+ D+ 


ie Ob Dye S © acl (J5DF > 0 


EY _ ( uk sD)0¥5D) (Kad Ex 
Dt (WED VPXKesDP Anh] | OF 





2 [ (WEsD)(Y3D) Kea DF |. 
T* UGaeap) (aD) kadytanh (DCI 


caontenrlckaoF- ooh EFC] 


CWED)CYDICKea0- jot 


[Keay=Ckes pod Ich. aye 
+ (eto cath 


Ber _ [emer aT er Peer) 
Dt = © (WD) CLD)Kesdo- ot 


+2 Choad CED) Ker DY = (Kes DY J 
(WMD) CY.P)CKesPy— 


3 ye 
eS 


4 (CaP ead)'taunh OX) | E cry 
C¥aD) (Kea D~ pt 


Bo Cy : esD)3D) (Ker DY 
p+? Fa ak ClbOT 
BD) (Ker DY - 4) = 1G ys ] 
Cu) DCD) (Ke3Dy~ D+ 


. (WEs0) 30) Kea) ea” 
” oes eit yot | 


e bag D Bem -Coth ty ct cy" 
JL @Eenwmnckorcant JUS pt 


_: la>(eD Lea GaPFI) ' 
J [ (WED) (YD) (Kes3D 


f (¥5D) (Kad) kanh (030)C%)) 7 | ? 
(¥2D)C Kea D)* p+ 


77 


(psy ph - Be. (BNF 
p+ Dt 


Cys 33D) Cd) (Kea 
pt IK 


| 


wE2D)(HD (Ke 3DP Awmhl (¥30)( mill Dt 


Sp GlnD(6D) [Kea DI-CKe30F J ce 
iL CWE2 PY ID)(KsDP | | 


= 


eye. ch eg pe. (chyt 
ofp 


pe. (EF _ [ (WE3P)CY3D)KeaD)_ (eo) 
Dt CWE2 DX LP (KD) Ah LOSE WIL pt 


+ oMiELOaP™ Cet) Ee 7 ) 


(WEDD) CKD jp 


CBee c& _ pe. 
ot Dt 


pri(eyr _ Ee )CkadDP~ Ckes0t] [e-eF 
pr (W MDCLD NK 307 p+ 


+) neta Seah EShC OY cor ch CF 


C¥20)(Ke3D)* 


cehyFc® BEF 
pt D+ 





Bt hyF cy. = [Ganxennt than GedIN) CF ef. ch ] 
CW) MD) CD) (Kes DY~ jy Dt 


+ (2 (Kad) Aanh (C¥5DX6)) H (ay er) 
CLD)( Kea DY” 


(ehyrch _ BA. (chy 
ia De 


rer: wonky fi stank Lou c6} et “( ct 


tanh? (3d) el 


oor ov = (EF ( |— tank TO, erro foe cy] 


roy. avr y teal Cb 
- 


ae aed ees 


otyint _ _ peor 
p+  =—s—i‘iéi 


REGION 1 
The term Pjq will be used for the power flow in region 1 
for the case (y,D)* > 0. 
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For the case (Y¥,D)* < 0, Yy, is imaginary in which case 
cy, Da = -(¥,D)*. The power flow, Pi, , for this case is 
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REGION 2 

For region 2, the power flow expressions are the same as 
for unshielded slotline except that tthe Fourier integral is 
replaced by a summation and the interval 2 is replaced by 
the interval b. 
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REGION 3 
For the case (Y; D)* oak is real and the power flow, 
P30 > 1s 
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For the case (¥3D)* < 0, ¥3 1s imaginary in which case 
(Y,D)* = -(Y3D)*. The power flow, Pap, > Loumbills Case 1s 
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